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Abstract. Direct imagery and long-slit, spatially re- 
solved echellograms of the high excitation planetary neb- 
ula NGC 1501 allowed us to study in detail the expansion 
velocity field, the physical conditions (electron tempera- 
ture, electron density, ionization) and the spatial distribu- 
tion of the nebular gas. 

An electron temperature of 11500 K and a turbulence 
of 18 km s~^ are derived by comparing the Ha and [OIII] 
emission line profiles, but large, small scale fluctuations of 
both these quantities are present in the ionized gas. 

The radial density distribution shows external peaks 
up to 1400 cm~^; they have steep outwards profiles and 
extended inwards tails probably originated by Rayleigh- 
Taylor instability and winds interaction. 

The complexity of the expanding motions indicates 
that the main part of NGC 1501 is a thin ellipsoid of mod- 
erate ellipticity, but the presence of a pair of large lobes 
along both the major and the intermediate axes and of a 
multitude of smaller bumps spread on the whole nebular 
surface, makes the general 3-D structure of NGC 1501 like 
a boiling, tetra-lobed shell. 

This peculiar morphology can be qualitatively ex- 
plained in terms of interaction of the slow nebular ma- 
terial with the intense and fast wind from the WC4/0VI 
central star. 




Fig. 1. Broad-band R image of NGC 1501 (exposure time 
600s, seeing 0.76") obtained with the 3.58m Italian Na- 
tional Telescope (TNG); North is up and East is on the 
left. At higher contrast a very faint, diffuse and roundish 
envelope appears, extending up to 34" from the central 
star. 



Key words: planetary nebulae: individual: NGC 1501 - 
ISM: kinematics and dynamics - ISM: spatial models 

1. Introduction 

The little studied, high excitation planetary nebula (PN) 
NGC 1501 (PNG144.5-h06.5, Acker et al., 1992) is "a 
very irregular and patchy elliptical disk, about 56"x48" 
in P.A.=98°. The periphery shows traces of a broken ring 
formation; the brightest portions are the edges at the ends 
of the minor axis. Relatively faint" (Curtis, 1918). 

At least a dozen of distance estimates of this nebula 
are contained in the literature. They span in the range 
0.9 Kpc (Amnuel et al., 1984; statistical distance based 
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on the surface radio-brightness) to 2.0 Kpc (Acker, 1978; 
individual distance based on the average extinction in the 
galactic disk) and roughly peak around 1.3 Kpc (the value 
we will adopt in the present paper). 

A kinematical study by Sabbadin and Hamzaoglu 
(1982 a) suggests that NGC 1501 is a prolate spheroid 
of moderate ellipticity. 

The spectral type of the exciting star is WC4/ OVI 
(AUer, 1976, Tylenda et al., 1993, Corny and Stasin- 
ska, 1995). The star has a high temperature, between 
logTz(Hcii) = 4.91 - 4.98 K (Sabbadin, 1986, Stanghellini 
et al., 1994) and logTmodd atmosphere = 5.13 K (Koesterke 
and Hamann, 1997a). It is loosing mass at a rate of 5.2 x 
10~^ Mo yr"-"^ (Koesterke and Hamann, 1997a) and ter- 
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Fig. 2. The same R frame of Fig. 1 after the appUcation 
of a soft Lucy-Richardson restoration (point spread func- 
tion=seeing, 10 iterations), showing in great detail the 
tangled Ha structure of NGC 1501. The haloes around 
stars projected onto the nebula are artifacts due to the 
deconvolution. 

minal wind velocities of 1800 or 3300 kms""'^ (Koesterke 
and Hamann, 1997a and Feibelman, 1998, respectively). 

It is one of few PNe nuclei showing nonradial g-mode 
pulsations (Bond et al., 1993, 1996, and CiarduUo and 
Bond, 1996). Following the "born-again scenario" pro- 
posed by Iben et al. (1983; see also Blocker, 1985 and Her- 
wig et al., 1999), these extremely hot, hydrogen deficient 
stars suffer a late thermal pulse (ejecting the hydrogen 
rich layers and exposing the naked C/0 core) after they 
reached the White Dwarf cooling sequence and are possi- 
ble precursors of the short-period GW Vir (PG 1159-035) 
pulsating White Dwarfs. Hamann (1997) and Koesterke 
and Hamann (1997b) suggest the evolutionary sequence 
[WC-late] ^ [WC-early] ^ [WC-PG1159] 
PG1159. 

Recently we started a long-term observing program 
aimed to investigate the physical conditions and the spa- 
tial structure of selected PNe through direct imagery and 
long-slit Echelle spectra. This article concerns NGC 1501 
and its plane is as follows: in Section 2 we present the 
observational material; Section 3 contains the expansion 
velocity field; the physical conditions (electron tempera- 
ture, turbulence and electron density) of the ionized gas 
are analysed in Section 4 and the tomographic maps in 
Section 5; in Section 6 we describe and discuss the result- 
ing spatial model; conclusions are drawn in Section 7. 

2. Observational material 



2.1. Imagery 

The best imagery of NGC 1501 available in the litera- 
ture goes back to Minkowski (1968), who obtained Ha 
and [OIII] interference filter plates with the Hale Tele- 
scope at Mount Palomar, revealing the intricate nebular 
structure (but see also Chu et al., 1987, and Manchado et 
al., 1996a). 

The superior spatial resolution of HST hasn't yet been 
exploited for NGC 1501; at the moment, only four broad- 
band, quite underexposed frames of the nebula (taken by 
Howard Bond on September 7, 1995, during a snapshot 
survey for companions of PNe nuclei) are contained in the 
HST archives; because of the low surface brightness of 
the emitting gas, they little add to Minkowski's Ha plate 
published in the frontispiece of the lAU Symposium N. 34 
proceedings. 

On September 1999, a broad-band R image of NGC 
1501 (exposure time 600s, seeing 0.76") was obtained with 
the Optical Imager Galileo (OIG) mounted on the 3.58m 
Italian National Telescope (TNG, Roque de los Mucha- 
chos. La Palma, Canary Islands) and equipped with a 
mosaic of two thinned, back-illuminated EEV42-80 CCDs 
with 2048x4096 pixels each (pixel size=13.5 microns; pixel 
scale in 2x2 binned mode=0.144 " pix'^). The CCD data 
array was processed in the usual way using standard IRAF 
packages. 

The appearance of NGC 1501 in the R band, shown 
in Figure 1, is almost entirely due to Ha emission; the 
contamination of [Nil] (AA6548 and 6584 A) and of Hel 
(A5876 A) amounts to 4% and 3%, respectively; other 
fainter lines contribute for less than 2% (see Kaler, 1976, 
Stanghellini et al., 1994, and the next sub-section). 

Our R image reveals the presence of a diffuse, roundish 
envelope framing the main nebular body and extending up 
to a distance of 34" from the central star (this halo is too 
faint to be reproduced in Fig. 1). 

The application of a soft Lucy-Richardson restoration 
(point spread function=seeing, 10 iterations; Figure 2) to 
the original R frame highlights the bubbly structure of 
NGC 1501, rich of filaments and condensations ("bear- 
ing a resemblance to the convolutions of the brain" is the 
suggestive description given by Francis G. Pease in 1917). 

2.2. High resolution spectra 

On December 1998, AA4500-8000 A spatially resolved, 
long-slit spectra of NGC 1501 {+ flat fields + Th-Ar cal- 
ibrations -|- comparison star spectra) were secured with 
the Echelle spectrograph (Sabbadin and Hamzaoglu, 1981, 
1982 b) attached to the Cassegrain focus of the 182cm Asi- 
ago telescope and equipped with a Thompson 1024x1024 
pixels CCD. 

We selected four position angles: 100° (apparent major 

axis, after Curtis, 1918), 10° (perpendicular to the appar- 
ent major axis), 55° and 145° (intermediate positions). 
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Table 1. Integrated line intensities 



A(A) 


ion 


I(obs.) 


I(corr.) 


4686 


Hell 


30 


33 


4861 


HI 


100 


100 


4959 


[OIII] 


430 


400 


5007 


[OIII] 


1300 


1200 


5876 


Hel 


16 


10 


6560 


Hell 


10: 


4: 


6563 


HI 


665 


285 


6584 


[Nil] 


24 


10 


7135 


[Arlll] 


22 


7 



All spectra were centred on the exciting star; we used the 
stellar continuum as a position marker and to correct the 
data for seeing and guiding uncertainties. The slit-width 
was 0.200 mm (2.5" on the sky), corresponding to a spec- 
tral resolution of 13.5 km s~^ (1.5 pixel). The calibration 
in wavelength and flux was performed in the standard way 
using the IRAF Echelle packages. 

Ha and A5007 A of [OIII], the dominant emissions in 
our spectra, present the same structure in every detail 
(but not in the line width: as normally observed, the first 
is broader than the second, because of the larger thermal 
motions). Besides Ha and A5007 A (and the correlated 
lines H/3 and A4959 A), only a few, faint emissions were 
detected, due to the low surface brightness of the emit- 
ting gas; they all mimic the Ha and [OIII] structure. The 
line intensities, integrated over the whole nebula and the 
entire velocity profile, are listed in Table 1; only relative 
fluxes are reported here, since the observational conditions 
(the sky was stable only at intervals) prevent the accurate 
calibration into absolute fluxes. 

The interstellar extinction was derived by comparing 
the observed Ha/H/3 intensity ratio with the dereddened 
value given by Brocklehurst (1971) for Te = 10^ K and 
Ne = 10'* cm-3; we obtain c(H/3) = 1.05 ± 0.10, in agree- 
ment with the values of 1.1 and 0.96 reported by Kaler 
(1976) and StanghcUini et al. (1994), respectively. 

The close resemblance of the Ha and [OIII] emission 
structure and the weakness (or the absence) of low excita- 
tion lines (in particular, [Nil] at AA6548 and 6584 A and 
[SII] at AA6717 and 6731 A) indicate that NGC 1501 is an 
optically thin, density bounded PN. 

3. The expansion velocity field 

The peak separation, measured at the centre of each neb- 
ular emission, gives the expansion velocities contained in 
Table 2, where ions are put in order of increasing ioniza- 
tion potential. For all ions (but not HI) the estimated error 
in the expansion velocity is essentially anti-correlated to 
the emission intensity; the large uncertainty in Ha comes 
from the line broadening due to thermal motions. 



Table 2. Expansion velocities 



A 


ion 


LP. 


2 Vexp 


(A) 




(eV) 


(km s"^) 


6563 


HI 


13.6 


80 ±3 


6584 


[Nil] 


14.5 


84 ±3 


5876 


Hcl 


24.6 


82 ±3 


7135 


[Arlll] 


27.6 


82 ±3 


5007 


[OIII] 


35.1 


81 ±2 


4686 


Hell 


54.4 


79 ±3 



Previous expansion velocity measurements in the neb- 
ula go back to Robinson et al. (1982, 2wcxp[0III]) = 78 
km s~^) and to Sabbadin and Hamzaoglu (1982a, 
2ve^p{Ra) = 2«exp([0///]) = 76 km s-*). 

The very low (if any) expansion gradient present in 
NGC 1501 confirms that stratification effects are negligible 
and suggests that the main nebular emission occurs in a 
narrow shell (note, in particular, the large Ha expansion 
velocity, as large as for the other lines). 

More precise informations can be obtained from the 
detailed analysis of the expansion velocity field at the four 
position angles. These are shown in Figure 3 for Ha; the 
velocity and intensity distributions derived for A5007 A of 
[OIII] coincide with the Ha ones; though the [Nil] line at 
A6584 A is too weak for an accurate study, its intensity 
and velocity trends appear very similar to those observed 
in Ha and [OIII]. 

The Ha radial velocity distributions of Fig. 3 present 
the classical bowed shape expected of an expanding shell, 
but they are so inhomogeneous and distorted that a simple 
model (like a triaxial ellipsoid) is decidedly inadequate to 
fit all the data. Note, in particular: 

- the comparable nebular extent at P.A.=10° (perpen- 
dicular to the apparent major axis) and at P.A.=55° (in- 
termediate direction); 

- the noticeable differences in both the expansion veloc- 
ity field and the nebular structure at P.A.=55° and 145° 
(being these two position angles symmetrically arranged 
with respect to the major axis, we would expect almost 
mirror structures for a simple rotational figure); 

- the similarity of both the velocity field and the in- 
tensity distribution at P.A.=10° and 145° (suggesting that 
one of the axes of symmetry is projected in P.A.~170°). 

Moreover, the imagery of the nebula (taken a few 
months after the spectroscopic material) stands out a fur- 
ther complication: the direction of the apparent major axis 
of NGC 1501 seems closer to 110° than to 100° (the P. A. 
we used), or to 120° and 98° (the P. A. given by Pease, 
1917, and Curtis, 1918, respectively). 

To make easier the interpretation of the velocity maps 
of Fig. 3 in terms of an acceptable model, we believe con- 
venient the introduction of a parting line at an apparent 
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Fig. 3. Ha position- velocity contour maps observed in NGC 1501 at P.A.=10° (almost perpendicular to the direction 
of the apparent major axis), P.A.=55°, P.A.=100° (close to the direction of the apparent major axis) and P.A.=145°. 
"Low" and "high latitude" nebular regions, as introduced in the text, are indicated. 



distance of about 15" from the central star, separating the 
"low latitude" ionized gas from the "high latitude" one. 

The "low latitude" velocities arc quite regular ellipti- 
cal arcs (i.e. projections of a triaxial ellipsoid), whose tilt 
is maximum at P.A.=10° and P.A.=145°, and minimum 
(but not null) at P.A.=100° (i.e. close to the apparent ma- 
jor axis). To be noticed that the un-tilted velocity ellipse 
occurs at P.A.~80° (corresponding to the hue of nodes; 
this agrees with the foregoing indication that one of the 
axes of symmetry is projected at P.A.~170°). The inten- 
sity distribution appears knotty and irregular at P.A.= 
55° and 100°, while at P.A.=10° and 145° it presents two 
symmetric, extended condensations (as normally observed 
in bipolar PNe). 

The "high latitude" velocity maps of Fig. 3 are charac- 
terized by an extreme variability: they are faint, extended 
and asymmetric at P.A.=10° and 145°, faint, extended 
and quite symmetric at P.A.=100°, knotty, bright, sym- 
metric and small at P.A.=55°. The general impression is 



that they represent hemispheric bubbles protruding from 
the ellipsoidal region of NGC 1501; in some directions (as 
in the East sector of P.A.=100°) multiple structures ap- 
pear, i.e. an external, smaller dome overlaps the internal, 
larger one. 

A further, remarkable feature of the "high latitude" 
zones in Fig. 3 is the internal, weak, diffuse emission 
present at all the four position angles, indicating that a 
low density ionized gas completely fills these nebular re- 
gions. 

In order to decipher the complex gas motion observed 
in NGC 1501, a detailed analysis of the nebular physical 
conditions will be performed in the next Section. 

4. Physical conditions 
4.1. Electron temperature (Te) 

In absence of diagnostic line intensity ratios (e.g. 
4363/5007 A of [QUI] and 5755/6584 A of [Nil]), the elec- 
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tron temperature of the ionized gas was derived by com- 
paring the Ha and [OIII] emission line profiles; the basic 
assumption is that, for a given Te, the thermal motion in 
H is four times larger that in O (the former element being 
sixteen times lighter that the second). 

CAVEAT: r/iis method can be rightly applied to PNe 
only if the and O^"*" layers do coincide ( as in NGC 
1501, which is an optically thin, high excitation PN). In 
many cases, the presence of large stratification effects in- 
validates the results, since and O"*""*" are emitted in 
separated layers expanding in different ways (due to the 
presence of a gradient in the velocity field). Similarly, the 
77+ and line profile comparison fails in medium and 
high excitation PNe. 

The full width at half maximum (W) of an emission 
line is the convolution of different components: 

^^total ^^instr. ~b ^^thermal ^^turbulence ^^gradient ("^) 

where the last (too frequently neglected) term takes 
into account the existence of a radial gradient in the ex- 
pansion velocity field. 

A further broadening factor concerns Ha: its seven 
fine structure components can be modeled as the sum 
of two equal Gaussians separated by 0.14 A (Dyson and 
Meaburn, 1971, and Dopita, 1972). 

Since in NGC 1501 the H+ and 0++ layers do coincide, 
we have: 

Te = f[W{Ha)l,,i - W{[OIII])l,J (2) 

A mean value of Te = 11500 ± 500 K is derived by 
analysing both the blue-shifted and red-shifted Ha and 
[OIII] profiles in the central region of the lines; the spread 
in Te, only partially due to instrumental -|- measurement 
uncertainties, testifies the existence of electron tempera- 
ture variations within the ionized gas, although no clear 
correlation was found with position and/or flux. 

The only previous Te determination in this nebula 
is reported by Stanghellini et al. (1994), who obtained 
Te=10700 K from the 4363/5007 A [OIII] intensity ra- 
tio. Moreover, some statistical works on a large sample 
of PNe (e.g. Cahn et al., 1992, and Phillips, 1998) adopt 
Te=15100 K, based on the strength of the Hell A4686 A 
line and Kaler's, 1986, calibration. 

Having obtained Te, the value of (W'^f^^t + VF^,.„^ ) 
can be inferred from relation (1); since in NGC 1501 the 
expansion velocity gradient is small (see Table 2), it can 
be neglected, allowing us to quantify turbulent motions in 
the ionized gas. We obtain for Wturbuience a mean value of 
18 km s~^ from both the Ha and [OIII] profiles, but small 
scale fluctuations (up to 10 km s""'^) are present. 

A comparison with previous data reported in the liter- 
ature for PNe gives poor results, because of our "caveat", 
of the scarce bibliography and of the different reduction 
methods. We recall the analysis performed by Gesicki et 



al. (1998) on seven PNe, indicating that the highest tur- 
bulent motions (15 km s~^) occur in M 3-15, a nebula 
powered by a WC 4-6 star. Although this result is weak- 
ened by the evidence that M 3-15 is ionization bounded (at 
least in some directions), we agree with the statement of 
these authors: "It is possible that nebulae with [WC] cen- 
tral stars have less regular velocity fields than the other 
PNe... The assumed very high turbulence may only be 
an approximation to a more complicated situation with 
strong velocity variations in radial direction" . 

4.. 2. Electron density (Ne) 

Also A^e diagnostic line ratios (e.g. 6717/ 6731 A of 
[SII], 3726/3729 A of [Oil], 4711/4740 A of [ArlV] and 
5517/5537 A of [ClIII]) are absent in our echellograms of 
NGC 1501. In order to derive the electron density, we will 
proceed in three steps: 

a) determination of the relative iVe radial trend in the 
"zero-velocity pixel column" (as defined by Sabbadin et 
al., 2000); 

b) transformation to absolute A''e values by means of 
a suitable calibrator; 

c) extension of the absolute Ne determination to the 
whole nebula's slice covered by the slit. 

Points a) and b) will be discussed here and point c) in 
the next Section, dedicated to tomography. 

Let's consider a long-slit, spatially resolved, high res- 
olution spectrum of a typical planetary nebula expand- 
ing at fexp- Following Sabbadin et al. (2000), the zero- 
velocity pixel column of an emission line represents the 
slice of nebula centred in the plane of the sky crossing the 
exciting star and having a depth along the line of sight 
t = rAv/vcxp, where r is the nebular radius and A?; is 
the pixel spectral resolution. As an example, in the veloc- 
ity maps shown in Fig. 3, the observed Ha zero-velocity 
column at each position angle is represented by a vertical 
strip centred at -|-36 km s~^ (systemic radial velocity of 
NGC 1501) and 9.0 km s"^ wide. 

At every position along the slit, the intensity in the 
zero-velocity column is proportional to NeNi; in the case 
of complete ionization and Te=constant, loc A''e^ (we im- 
plicitly assume the constancy of the local filling factor, e/, 
within the slice of nebula identified by the zero-velocity 
pixel column; e/ is the fraction of the local nebular vol- 
ume which is actually filled by matter with density A^e; 
the local nebular volume is given by: (pixel area) xf, where 
t = rAv/vexp). 

In short: the zero- velocity pixel column, isolating a 
slice of nebula unaffected by the expansion velocity field, 
establishes a direct link between the intensity profile and 
the ionic and electron density distributions, thus allowing 
us the detailed analysis of the radial gas structure (and 
ionization) in the expanding nebula. 

The Ha zero-velocity column distributions, corrected 
for contamination by the adjacent spectral columns and 
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Fig. 4. Relative (left ordinate scale) and absolute (right 
ordinate scale) radial density trends obtained from the 
Hq zero- velocity pixel columns. Note the steep outwards 
profile and the broad inwards tail present at all the four 
position angles (less evident in P. A. =55°). 

for seeing and guiding uncertainties (for details, see Sab- 
badin et al., 2000), are shown in Figure 4; they are nor- 
malized to the strongest intensity (i.e. the peak in the N-E 
sector of P.A.=55°). 

The radial density profile of NGC 1501 results quite 
complex; the main features, common at all the four posi- 
tion angles, are the following: 

- the outermost parts correspond to the faint, roundish 
halo detected in our R image (representing the vestiges of 
the photospheric material ejected in the AGB phase); 

- the density peaks are located in the external nebu- 
lar regions (close to the edge of the bright disk visible in 
Figs. 1 and 2); their distribution confirms that the main 
component of NGC 1501 is a shell of moderate thickness; 

- the density peaks show steep outwards profiles; 

- inwards tails of different shapes are present at all the 
four position angles; their extent seems anti-correlated to 
the height of the density peak: small at P.A.=55°, interme- 
diate at P.A.=10°, broad at P.A.=145° (note the detached 
structure in the N-W sector) and at P.A.=100° (note the 
jagged profile). 



Very similar radial density trends arc obtained by 
analysing the [OIII] zero-velocity pixel columns. 

In order to scale to absolute values the relative elec- 
tron density profiles shown in Fig. 4, we will utilize the 
direct intensity- A'^e (cm~'^) relation given by the Ha zero- 
velocity pixel column of NGC 40 (observed in the same 
nights and with the same instrumental setup used for 
NGC 1501). 

We recall that the [SII] red doublet is quite strong in 

the low excitation PN NGC 40 (AUer and Epps, 1976, and 
Clegg et al., 1983), thus an accurate determination of Ne 
can be obtained from the diagnostic ratio 6717/6731 A 
(Sabbadin et al., 2000). 

From the general expression: 



AnD^FiHa) = hiy az,2N{H+)NeeATxr'^dr (3) 
Jo 

we obtain the following relation between the surface 
brightnesses of the zero-velocity pixel columns of NGC 
1501 and NGC 40: 



iVei 



1,2 



Jlt'cxpl?'2-P2£, 
-f2l'exp2ri 1 



0.5 



TC2 

Tei 



(4) 



where the suffixes (1) and (2) refer to NGC 1501 and 
to NGC 40, respectively, / is the Ha intensity in a pixel of 
the zero- velocity column (corrected for interstellar extinc- 
tion), r" is the apparent nebular radius, D is the distance 
and ei is the local filling factor (as previously defined). For 
both nebulae we assumed Ne = 1.1A''(H+). 

The calibration through relation (4) gives for NGC 
1501 the absolute radial density profiles shown in Fig. 
4, right ordinate scale. The following parameters (some 
taken from the literature) were adopted: 

NGC 40: c(H/?)=0.60, Te=8000 K, Uexp=25 km s-\ 



r"=20", 15=1.1 Kpc. 

NGC 1501: c(H/3)=1.05, Te=11500 K, t;exp=40 km 
s-i, r"=28", Kpc. 

Moreover, we assumed e,(NGC 40) = e((NGC 1501) = 
1; in general, this appears to be a reasonable condition; 
some particular cases will be discussed later-on. 

The estimated inaccuracy in the A^e scaling factor is 
±20%, mainly due to our poor knowledge of the PNe dis- 
tances; fortunately, in relation (4) the electron density has 
a low dependance on the distance. 

We wish to point out that the foregoing method [i.e. 
the use of a PN (NGC 40) as a density calibrator for an 
other PN (NGC 1501)] represents a procedural escamo- 
tage : we were forced to adopt it since the night sky during 
the observational run was stable only at intervals. A more 
elegant analysis, based on the "absolute" surface bright- 
ness [i.e. formula (3)] will be presented and discussed in a 
future paper dedicated to the double envelope PN NGC 
2022. 
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Fig. 5. Spatial reconstruction of the Ne structure of NGC 1501 derived at P.A.=10° (almost perpendicular to the 
direction of the apparent major axis), P.A.=55°, P.A.=100° (close to the direction of the apparent major axis) and 
P.A.=145°. The minimum density shown is A'^e=200 cm~"^, while the maximum density (A^c=1380 cm~'^) is reached 
by the approaching gas located in P.A.=145° (N-W sector), at an apparent distance of 13" from the central star. Being 
NGC 1501 an optically thin, high excitation PN, the A''(H+) and A''(0++) tomographic maps coincide with the Ne 
ones (but for a scaling factor of 1.1 and 3.67 x 10^, respectively). 



5. Tomography 



To derive the de-projected expansion velocity field and the 

real gas distribution within the slice of nebula covered by 
the spectroscopic slit, we will follow the original method 



described by Sabbadin (1984) and Sabbadin et al. (1985), 
and recently applied by Sabbadin et al. (2000) to the low 
excitation PN NGC 40. 



The basic consideration is that the position, radial 
thickness and density of each emitting region in an ex- 
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tended, expanding object can be obtained from the veloc- 
ity, FWHM and flux, respectively. 

In the case of NGC 1501, the limited number of emis- 
sions detected and the total absence of stratification effects 
across the nebula force us to use a linear position-speed 
relation. Moreover, from the angular extent of the ellipses 
fitting the "low latitude" nebular regions in Fig. 3, we will 
adopt ri — 20" , where ri is the nebular radius in the ra- 
dial direction at the apparent position of the central star 
(for details, see Sabbadin, 1984). Finally, an electron tem- 
perature of 11500 K and a turbulence of 18 km s~^ (both 
constant over the nebula) are assumed. 

The adoption of other (reasonable) sets of parameters 
doesn't substantially modify the results here obtained. 

Our tomographic analysis of NGC 1501 mainly con- 
cerns the electron density; the A^e structme of the neb- 
ula at the four position angles is shown in Figure 5. Ac- 
cording to the discussion above, the A''(H+) tomographic 
maps coincide with the A''e ones (but for a scaling factor of 
1.1). The same situation occurs for A^(0+"'"), but in this 
case the scaling factor is 3.67 x 10^; in fact, the A5007 
A intensity reported in Table 1 gives A''(0++)/A(H+) = 
3.0(±0.4) X lO"'^ (see AUer and Czyzak, 1983 and Aller 
and Keyes, 1987). Moreover, the ionization structure of 
NGC 1501 suggests that most oxygen is doubly ion- 
ized, so that N{0)/N{B.) ~ N{0++)/N{B.+ ); to be no- 
ticed that Stanghellini et al. (1994) report for this neb- 
ula N{0++)/N{B+) = 3.72 X 10"^ and Ar(0)/A'(H) = 
5.7 X 10"^. 

The only previous Ne determination in NGC 1501 

dates back to Aller and Epps (1976), who observed at 
low spectral resolution a small region centred 20" from the 
star (in P.A.=125°); they derived /(6717)/7(6731) = 0.78, 
corresponding to A''e=1200 cm'^ (for Te=11500 K). Their 
direction is intermediate between our P.A.=100° and P.A. 
=145"; at a distance of 20" from the central star we obtain 
the following density peaks: 

iVe(blue shifted)=900 cm'^ and Are(red shifted)=1000 
cm-3 in P.A.=100''; 

A'^e(blue shiftcd)=700 cm'^ and A''e(red shifted)=900 
cm-3 in P.A.=145°. 

Unfortunately, a more detailed comparison appears 
hazardous, due to the differences in the spectroscopic res- 
olution and reduction procedure, and, in particular, to 
the presence of small scale density fluctuations within the 
nebula. 

The de-projected expansion velocities of NGC 1501 
(directly obtainable from Fig. 5, given the linear position- 
speed relation used) span in the range 38(±2) to 55(±2) 
km s~^. The slowest motions occur in the densest re- 
gions at P.A. =10" and P.A. =145°; the combination (high 
density + low expansion velocity) here observed suggests 
that the minor axis of the central ellipsoid is projected 
at P.A.~170° (in agreement with the indications given in 
Section 3). The largest expansion velocities [55(±2) km 
s^^] correspond to the high latitude, untilted, hemispheric 



bubbles at P.A. =100°; in this case, the combination (low 
density -|- high velocity) suggests that we are (almost) 
observing the projection of the major axis of the central 
figure. Finally, from geometrical considerations, the pro- 
jection of the intermediate axis can be put at P.A.~30°. 

Figure 5 summarizes most of the observational results 
given in the previous sections, confirming the composite 
structure of the nebula: NGC 1501 is an ellipsoid of mod- 
erate ellipticity, deformed by a pair of large lobes along 
both the major and intermediate axes and by a number of 
minor bumps scattered on the whole nebular surface (in a 
few cases, see for instance the big "ear" in the N-W sector 
of P.A. =145", the dimensions of these "minor bumps" ap- 
pear comparable to those of the lobes related to the axes 
of the central ellipsoid). 

In Fig. 5, the absence of the broad, inwards tail at low 
latitude is probably due to instrumental limitations: intu- 
itively, our echellograms of an extended PN like NGC 1501 
have a "spatial" resolution along the slit which is better 
than the "spectral" resolution along the dispersion. Due to 
projection effects, the "low latitude" zones (dominated by 
the expansion velocity field) have a "spectral" resolution, 
the "zero- velocity pixel column" has a "spatial" resolu- 
tion (being unaffected by the expansion velocity field) and 
the "high latitude" zones an hybrid resolution. Clearly, 
this gradual variation in resolution conditions the tomo- 
graphic reconstruction. In other words: we cannot exclude 
the presence of a low density, inwards tail also at low lati- 
tude. Only deep observations at much higher spectral res- 
olution could solve the question. 

Having said this, we believe that Fig. 5 adequately re- 
produces the true matter distribution in NGC 1501, and 
that the faint, inwards emissions visible at high latitude 
do represent the trace of the original ellipsoidal structure 
(note, in particular, the sharp radial profile in the high lat- 
itude zone at P.A. =55° and the detached structure in the 
N-W sector of P. A. =145°, fitting the ellipsoid projection). 
This implies that some accelerating agent partly swept-up 
the lower density regions of the triaxial ellipsoid, causing 
both the extended, hemispheric bubbles and the broad, 
inwards density tails. The final result of this acceleration 
is constituted by the spatial structure illustrated in the 
next Section. 



6. Spatial model 

An opaque sketch of the resulting 3-D model for NGC 1501 
is given in Figure 6, superimposed to the R frame of the 
nebula. The outermost contour of the model is enhanced; 
moreover. Fig. 6 contains: 

- the directions (A, B and C) of the semi-axes a, b and 
c of the central ellipsoid; 

- the large lobes associated to the major and interme- 
diate axes; 
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Fig. 6. Opaque sketch of the resulting spatial structure of 
NGC 1501 (superimposed to the red frame of the nebula). 
The model contains: - the directions (A, B and C) of the 
semi-axes (a, b and c) of the central ellipsoid; - the large 
lobes associated to the major and intermediate axes; - 
some secondary lobes visible at (or near) the nebular edge; 
- the extent of the faint, roundish envelope detected in the 
R image. Same orientation and scale as Figs. 1 and 2. 



- some "secondary" bubbles identified in the spectra 
and/or in the imagery; duo to projection effects, only those 
at (or near) the nebular edge are detectable: 

- the extent of the faint, roundish, external envelope 
visible in the Ha image. 

Although rare, the pecidiar morphology of NGC 1501 
is not unique amongst PNe. A quick look at the main im- 
agery catalogues allowed us to identify three more candi- 
dates: A 72 and NGC 7094 are faint, high excitation PNe 
presenting a complex filamentary structure (Manchado et 
al., 1996a); IC 4642 is quite bright, at high excitation, 
decidedly tetra-lobed in both the Ha and [OHI] images 
published by Schwarz et al. (1992). 

Moreover, Manchado et al. (1996b) introduced the 
morphological class of quadrupolar PNe, containing five 
compact objects (M 2-46, K 3-24, M 1-75, M 3-28 and M 
4-14); these nebulae present two pairs of lobes, each pair 
symmetric with respect to a different axis. A sixth candi- 
date, NGC 6881, was added by Guerrero and Manchado 
(1998). Following these authors, a quadrupolar PN can 
be formed by precession of the rotation axis of the central 
AGB star, possibly in the presence of a binary companion, 
associated with multiple shell ejection. 

Finally, the general spatial structure derived for NGC 
1501 (i.e. two pairs of bipolar lobes on axes having dif- 
ferent orientations and intersecting at the central star) 



presents noticeable analogies with the model suggested by 
Bahck and Preston (1987, their Fig. 4) for NGC 6543. 

In the case of NGC 1501, the modest departures from 
the spherical symmetry can be explained in terms of small 
inhomogeneities occurred during the super-wind ejection 
(for instance, due to stellar rotation), later enhanced by 
photoionization and winds interaction (see Dwarkadas and 
Balick, 1998 and Garcia-Segura et al., 1999). The intro- 
duction of a close companion (as suggested for NGC 6543 
by Balick and Preston, 1987, and Miranda and Solf, 1992), 
or even of a binary companion (as proposed by Manchado 
et al., 1996b, for quadrupolar PNe) appears unnecessary 
for NGC 1501. 

Clearly, the tetra-lobed shape is a simplification of the 
true spatial structure of this nebula; a close inspection 
to the observational data indicates that each macro-lobe 
is constituted of a heap of small components, that some 
morphological diff'erences exist amongst the lobes associ- 
ated to the major axis of the ellipsoid and those connected 
with the intermediate axis, and that these lobes are only 
roughly aligned with the axes of the central figure. An ex- 
tended spectroscopic coverage of the nebula is in progress, 
in order to obtain the detailed spatial structure of the ion- 
ized gas. 

The most intriguiging characteristic of the matter dis- 
tribution in NGC 1501 is the presence of an inwards tail 
in the radial density trend; this tail, particularly evi- 
dent in the directions of the lobes, can be the result of 
hydrodynamic processes in the nebular shell. Following 
Capriotti (1973; see also Breitschwerdt and Kahn, 1990, 
Kahn and Breitschwerdt, 1990, and Garcia-Segura et al., 
1999), in the first evolutionary phases, when the PN is still 
ionization bounded, Ray leigh- Taylor instabilities occur at 
the ionization front, forming a series of knots, condensa- 
tions and radially arranged fingers (see Dyson, 1974, and 
Bertoldi and McKee, 1990) which expand slower than the 
ionization front. 

Similar instabilities are produced also by the interac- 
tion of the fast stellar wind with the low velocity nebular 
material (Vishniac, 1994, Garcia-Segura and Mac Low, 
1995, and Dwarkadas and Balick, 1998); in this case the 
optical thickness of the nebula is unimportant. If con- 
firmed, the extremely large value of M (logM = —6.28 
M0 yr-i) derived for the WC4/0VI nucleus of NGC 1501 
by Koesterke and Hamann (1997a) using the standard 
atmosphere model for Wolf-Rayet stars, would indicate 
winds interaction as the main responsible of both the ra- 
dial density distribution and the large expansion velocity 
observed in this nebula. 

To be noticed that, besides the dynamical effects on 
the nebular gas, an intense and lasting mass- loss of hy- 
drogen depleted and He. C and O enriched photospheric 
material would produce enhanced ionic and chemical com- 
position gradients across the nebula (as recently observed 
by Sabbadin et al., 2000 in NGC 40, a low excitation PN 
powered by a WC8 star). 
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Moreover: if winds interaction and/or Rayleigh- Taylor 
instabilities are the sources of the inwards tails detected 
in the density distribution of NGC 1501, these tails 
arc essentially constituted of knots and condensations; if 
they survived ionization and/or heating by conduction, 
we expect ei(tail) < ei{ma.m shell), where ei is the lo- 
cal filling factor, as previously defined. In other words: 
the true electron densities in the blobs and condensa- 
tions forming the inwards tail are larger by the factor 
[e;(main shell)/ej(tail)]'^-^ than the values shown in Figs. 
4 and 5, where we assumed ej(tail) = ej(main shell) = 1. 

Unfortunately, our observational material is inade- 
quate to transform the previous speculations into quan- 
titative results to be compared with the theoretical pre- 
dictions. Detailed, deep studies at higher spatial and spec- 
tral resolution of this (overlooked so far) PN are needed 
to answer the stimulating questions here excited. 

7. Conclusions 

The radial density distribution of the ionized gas in the 
high excitation PN NGC 1501 presents density peaks 
up to 1400 cm~^; they have steep outwards profiles and 
extended inwards tails probably originated by Rayleigh- 
Taylor instability and winds interaction. 

By comparing the Ha and [OIII] emission line profiles 
we derive an electron temperature of 11500 K and a tur- 
bulence of 18 km s~^, but large, small scale fluctuations 
of both these quantities are present in the ionized gas. 

The complex expansion velocity field observed in this 
nebula is fitted by an ellipsoid of moderate ellipticity, de- 
formed by a pair of large lobes along both the major and 
intermediate axes and by a multitude of minor bumps 
scattered on the whole nebular surface. 

The peculiar morphology of NGC 1501 can be qualita- 
tively explained in terms of interaction of the nebular gas 
with the intense and fast wind from the WC4/0VI central 
star. 

The exquisitely observational nature of this paper ex- 
empts us for hunting sophisticate evolutionary models; 
since our race is finished, we strecth the baton to someone 
else (hoping to see a "modellist" ready to catch it and to 
put on a winning burst of speed). 
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